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CONTRIBUTIONS FROM THE JEFFERSON PHYSICAL 
LABORATORY, HARVARD COLLEGE. 

ON THE CONDITIONS TO BE SATISFIED IF THE SUMS 
OF THE CORRESPONDING MEMBERS OF TWO PAIRS 
OF ORTHOGONAL FUNCTIONS OF TWO VARIABLES 
ARE TO BE THEMSELVES ORTHOGONAL. 

By B. Osgood Peiecb. 

Presented March 14, 1906. Received April 12, 1906. 

If fa (x, y), fa (x, y) are the potential functions due to two colum- 
nar distributions of matter the lines of which are perpendicular to 
the xy plane, and if ^i (x, y), fa (x,y) are conjugate to fa and fa, 
respectively, the families of curves obtained by equating ^i and fa to 
parameters, are lines of force of the two distributions. Moreover, 
fa + fa is the potential function due to a combination of the two 
distributions, and the function fa + fa equated to a parameter gives 
the corresponding lines of force. The fact that if (fa, fa) are any 
pair of conjugate functions and (fa, fa) any other such pair, the 
functions (a fa + b fa, a fa + bfa) are also conjugate — with similar 
facts for other classes of functions — lies at the foundation of the 
graphical methods so successfully used by Maxwell 1 and by others in 
drawing equipotential lines, and lines of force or flow, due to combi- 
nations of simple elements. If (fa, fa) are merely a pair of orthogonal 
functions and (fa, fa) another such pair, it is generally not true that 
(fa + fa, fa + fa) are an orthogonal pair : thus (x, y), (x* + y\ y/x) 
are pairs of orthogonal functions, but x + x 2 + y\ y + y /x are not 
orthogonal. 

In certain classes of physical problems one encounters potential 
functions which are not themselves harmonic and the lines of which 
are not possible lines of any harmonic function, and it is often de- 

1 Maxwell, Treatise on Electricity and Magnetism, Vol. I, Ch. VII. Minchin, 
Uniplanar Kinematics, § 112. See also P. W. Bridgman, The electrostatic field 
surrounding two special columnar elements, These Proceedings, 41, 28. 



150 PROCEEDINGS OF THE AMERICAN ACADEMY. 

sirable in cases where the analytical processes become too complex, 
to determine graphically the forms of lines of force or flow due to 
a combination of two simple elements. This note discusses briefly 
the conditions under which the ordinary method of procedure is 
possible. 

Let (a, /?) and (a, /*) be two pairs of orthogonal functions of the 
two variables (x, y), so that 

dx dx dy dy ' 

dX du. dX du. 

— ' — H • — = ()■ (2) 

dx dx dy dy ' v ' 

then if (a + A, /? + ft) are to form an orthogonal pair, the equation 

(da d±\m S ± \ (da dX\(dJ_ ,M_ (s) 
\dx ^ dx) \dx + dx) + \dy + dy) \dy dy J ~ V w 

must be identically satisfied. Since (1) and (2) are true, (3) takes 
the form 

/dk . dj$_ dX . d/3\ (da .da da . <V\ _ 

\dx dx dy dy J \dx dx dy dy) ~ *• ' 

If h a , h p , h x , h^ represent the values of the gradients of a, /?, A, /*, 
and if the angle at any point between the directions in which X and /? 
increase most rapidly be denoted by [a, /J], (4) becomes 

h K ■ hp • cos [X, /3] + h a ■ h^ ■ cos [a, p.] = 0. (5) 

Whatever the sequence of the directions of the gradient vectors 
might be, the two angles which appear in (5) would be either equal 
or supplementary, and their cosines would be equal in absolute value, 
but the gradients themselves are intrinsically positive and the sequences 
must therefore be such that 

**/*„ = Klh • ( 6 ) 

Suppose that in the case of two given pairs of orthogonal functions 
(a, /3,) (a, h), the necessary condition (6) is satisfied, and that the 
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value of the gradient ratio, hjh^, obtained from the given values of 
a and /3, is the function Q of x and y ; then 

and if, for (do.y/(dy) 2 in this equation, we substitute the value 

mmm 

obtained from (1), it appears, since the gradient of a real function 
cannot vanish, that 



da 
ox 



(I)- 

If for {daflidxf in (7), we substitute its equivalent 



(10) 

derived from (1), we shall learn that 

Either the upper signs or the lower signs must be used in (9) and (11). 

If now we treat the equation 

V = « 2 V (12) 

in a similar way we shall obtain the equations 



^= ±Q K<§)> (13) 



dx 



ax 

dy 



-"»(£)• <I4) 
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and, so far as the relation (12) is concerned, we may use either the 
upper signs or the lower signs, but if (4) is to be identically satisfied, 
the same sign must be used in (9) and (13) and the sign opposite to 
this in (11) and (14). Equation (6), then, together with the proper 
choice of sequence of directions for the gradient vectors which corre- 
sponds to the convention with regard to signs just made, will ensure 
the orthogonality of a + X, fi + /*. For practical purposes, however, 
it is well to approach the problem from another side. 

If (a, /?) and (X, p) are given pairs of orthogonal functions, and if 
we denote the given scalar point functions obtained by dividing d/3/dx 
by da/dy, and by dividing dp/dx by dk/dy, by f and t\, the equations 
(1) and (2) can be written in the forms 

dx\ dy oyl dx 

= V, (16) 



and 

djildX 

dy ~ dyj dx 



djLJd\_ 

dxj dy ~ 



or 



and 



dx~*° dy' dy~~ 4 dx' ^ } 



dft. 8k dfL 8k . . 

dx = 7 >'dy' Ty-'^'Tx- (18) 



If the values of the derivatives of 1$ and /* given in (17) and (18) 
be substituted in (4) this equation becomes 

«-«(b-!J-b-is)-* (19) 

and if (a + A, fi + /x) are to be orthogonal, a and X must be such as 
to satisfy it. If X were expressible as a function of o, and /a as a 
function of /?, the second factor would vanish, but this case is of no 
practical interest and (19) demands in general that £ and rj shall be 
identical, so that 
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and 



d$_ Ba ?/?__, .^ ,.-s 

Bx~^ By' By~ 4 Bx' y J 

dj._ t .B\ dp _ . . _ 3X , 20 v 

Bx~^ By' By c ~ Bx' K } 



If in these equations the arbitrary function £ is made equal to unity, 
the conditions degenerate into the familiar definitions of any two 
pairs of conjugate functions. 

In order that a single function (/3) may exist the partial derivatives 
of which with respect to x and y shall be equal, respectively, to 

t Ba , Ba 

V Ty^~ V Yx' 

it is necessary and it is sufficient that a and £ should satisfy the 
condition 

B f t 8a\ , B /, Ba\ „ 
By{^B-y) + rx{^Bx)=°' ^ 

or 

In order that /* may exist, £ and X must satisfy the equation 

S-B + g-g + f-W-' W 

If loge £ be represented by ar, the last two equations take the 
forms 

dm Ba Bsr . ^ a , »« / \ _ n / \ 

Bx Bx By By ^ ' ' *■ ' 

Sir _ 8A. Bm _ 9 A T72/\\ A 

Bx Bx By By ^ ' ' ^ ' 

and, if each of these be differentiated with respect to x and with 
respect to y, m may be eliminated from the resulting equations and a 
necessary condition for a and a obtained, which may be stated in the 
form of the determinantal equation — 
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da 
dx 


da 
dy 











V(a) 
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dx 


dx 

dy 
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dx ' dy 
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d 2 X 

dx* 
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dx 
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d 2 a 
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d* a 
df 
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a a A 

dx ' dy 


d 2 X 
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dX 
8x 


dX 


riM 



0. (26) 



If a and A happen to be harmonic, the elements of the last column 
vanish and the equation is satisfied, as it should be. 

It is possible to factor the determinant, after it has been reduced, 
and if 



T —^a L ,dX_d^X 4 da 



dx 



d 2 X da 



dx 2 dx dx 2 dx dx ' dy dy dx ' dy dy' 



M = 



d 2 a dX 



3 2 A da 



■ aT + 5T3 * 37 ~ • 



dX d 2 X da 



dx ' dy dx dx ' dy dx dy 2 dy dy 2 dy ' 



dX 

*y 



(27) 



-hi**)-* 

the condition of (26) demands either that a and X satisfy the 
equation 

V»(a) 



da 
dx 


da 
dy 


8X 

dx 


dX 
dy 



M 



V(A) 

N 



0, 



(28) 



PEIBCE. — THE GRAPHICAL SUPERPOSITION OF LINES OF FORCE. 155 



or else that o and X have the same level curves: this last case, as 
being uninteresting, may be left out of account. 

Sometimes (28) is more convenient than the unexpanded form of 
the same condition which follows immediately if we solve (24) and 
(25) for dar/dx and da/dy, 



(29) 



9n7 


_S"" w - 


JX 

dy 


V 2 (a) 


dx 

3s7 


da dX 
dx dy 


da 
dy 

da 

dx 


dX 

dx 

V°(X) 


dy 


da dX 
dx dy 


da 
dy 


ax 

dx 



(30) 



and equate the derivative with respect to y of the second member of 
the first equation to the derivative with respect to x of the second 
member of the other. This process yields the relation, 



_a_ 

dy 



da d\ da _ dX 
dx dy dy dx 



d_ 

dx 



■dX 

dx 



•V»(«). 



da 

dx 



V*(X) 



da dX 9o i dX 
dx dy 9y dx 



(31) 



and it is possible to check the fact that (28) and (31) are equivalent 
by a straightforward but somewhat laborious comparison of the two. 

If a and X satisfy (31), a function £ exists which satisfies (22) and 
(23), functions /? and /i exist which satisfy (17) and (20), and (a, /3), 
(X /*), (a + X, /} + ft) are orthogonal pairs of functions. 

If, for instance, both a and X represent values in the xy plane of 
harmonic space functions ( V, W) the level surfaces of which are sur- 
faces of revolution about the x axis, so that 



d_ 
dy 



( dV\ d*V n 



with a similar equation for W, 



"w=-^. "»- 



1 ,8X 



(32) 



(33) 
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equation (31) is satisfied and 

dsr . da 1 , . 

^ = 0, ly~ = y> l = ey - (34) 

In this case, if we put c = 1, /? and n are the Stokes functions corre- 
sponding to a and A. If the level surfaces of the harmonic space 
functions, V and W, are surfaces of revolution about two different 
straight lines in the xy plane, the functions a and X which represent 
the values of V and W in this plane do not in general satisfy (31). 

Graphical superposition of the lines of force in the xy plane due to 
an infinitely long, homogeneous cylinder of revolution parallel to the 
axis, and to a homogeneous sphere with centre in the plane, will not in 
general yield the lines of force in the xy plane due to a combination of 
the two masses. 

If a and X are harmonic, any linear function (but no other than 
a linear function) of o is harmonic, and any two linear functions of 
o and A. satisfy (31). There generally exist, however, non-linear 
functions of a. and X which, although they are not harmonic, satisfy the 
condition. The functions (x 2 — y 2 ), (x 2 + y 2 )", the second of which is 
not harmonic, obey (31), as do the harmonic pair (x 2 — y 2 ), log (x 2 + y 2 ). 

As a simple example of the fact that a harmonic function and a 
function which is not even isothermal may satisfy the condition (31), 
we may consider (2y* — x 2 ) and (y 2 — x 2 ). 

- a x 2 , which are solu- 

= 0, (35) 

"•» y "i> 

evidently satisfy the equation (31). 
If o and X are any two solutions of the equation 

d 2 V d 2 V dV 

i£ + W + SV>Z = 0. (36) 

where f(x) is any given function of x, the condition (31) is satisfied 
and or =/(#). 



The non-isothermal functions 
tions of the equation 
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If a and A satisfy the equation 

w+W + a ^ + b "¥ = 0, (37) 

st is of the form a x + by + c. 

In general a and X must both be solutions of an equation of the 
form 

where P and Q are any functions of x and y such that dPJdy = dQ/dx. 

The question whether if (a, y8) and (A, /x) are orthogonal pairs and 
(o + A, f$ + /i) is not an orthogonal pair, it is possible to find a func- 
tion (i?) of /?, and a function (M) of //. such that (a + A, £ + M) 
shall form an orthogonal pair, has already been answered ; for a and X 
must satisfy (31) in any case, and if they do this, or may be determined 
from (29) and (30) and B and if from (17) and (20). 

The Jefferson Laboratory, 

Harvard College, Cambridge, Mass. 



